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Abstract 

Observations have recently indicated that the Universe at the present 
stage is in an accelerating expansion, a process that has great implica- 
tions. We evaluate the spectrum of relic gravitational waves in the current 
accelerating Universe and find that there are new features appearing in 
the resulting spectrum as compared to the decelerating models. In the 
low frequency range the peak of spectrum is now located at a frequency 
i/E — {%^y^^i^H^ where is the Hubble frequency, and there appears 
a new segment of spectrum between ue and uh- In all other intervals 
of frequencies > the spectral amplitude acquires an extra factor 
due to the current acceleration, otherwise the shape of spectrum is simi- 
lar to that in the decelerating models. The recent WMAP result of CMB 
anisotropics is used to normalize the amplitude for gravitational waves. 
The slope of the power spectrum depends sensitively on the scale factor 
a(r) oc |r|^+^ during the inflationary stage with /3 = — 2 for the exact de 
Sitter space. With the increasing of the resulting spectrum is tilted 
to be flatter with more power on high frequencies, and the sensitivity of 
the second science run of the LIGO detectors puts a restriction on the 
parameter (5 < —1.8. We also give a numerical solution which confirms 
these features. 
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1. Introduction. 



The inflationary expansion of the early Universe can create a stochastic background of relic gravita- 
tional waves, which is important in cosmology and has been extensively studied in the past [1] [2] [3] 
[4] [5]. The spectrum of relic gravitational waves, as is to be observed today, depends not only on the 
details of the early stage of inflationary expansion, but also on the expansion behavior of the subsequent 
stages, including the current one. The calculations of spectrum so far [6] [7] [8] [9] [10] [11] have been 
done for the case that the current stage is in a decelerating expansion. The resulting spectrum has 
been put among the candidate list of sources for the gravitational wave detectors either in operation 
[12] [13], or under construction [14] [15]. The astronomical observations on SN la [16] [17] indicate that 
the Universe is currently under accelerating expansion, which may be driven by the cosmic dark energy 
{Qa ~ 0.7) plus the dark matter ~ 0.3) [18]. This is further supported by the recent WMAP results 
on CMB anisotropics [19] with Q = 1. By the wave equation, the evolution of relic gravitational waves 
depends on the expanding spacetime background, and the wave amplitudes depend on whether the wave 
lengths are inside or outside the Hubble radius. The Universe under accelerating expansion has a Hubble 
radius as a function of time that differs from that in the conventional models of decelerating expansion. 
Therefore, one can expect that, if an earlier matter-dominated stage of decelerating expansion is followed 
by the current accelerating expansion, the outcome for the spectrum of relic gravitational waves will be 
altered. 

In this paper we study the impact of the current accelerating expansion on the relic gravitational 
waves. We will first sketch, as a setup, the well-known formulations of the gravitational waves in an 
expanding spacetime [20], and give the explicit scale factor a(T) for a sequence of successive expanding 
epochs, including the current epoch of accelerating expansion. We then evaluate the power spectrum 
of relic gravitational waves and find that the current accelerating expansion does change the spectrum, 
including its shape and amplitude. Throughout the paper we shall work with units with c — h — 1, 
otherwise it will be pointed out. We also use notations similar to that of Grishchuk [8] for convenience 
for comparison. 



Incorporating the perturbations to the spatially flat Robertson- Walker space-time, the metric is 



where where r is the conformal time, the perturbations of spacetime hij is a 3 x 3 symmetric matrix 
containing generally the scalar, vector, and tensor parts. The gravitational wave field is the tensorial 
portion of hij, which is transverse-traceless 



2. The Gravitational Wave Equation. 



ds^ = a^{T)[dT'^ — {Sij + hij)dx^dx^], 



(1) 



dih'^ = 0, 5'^ hij = 0. 



(2) 



We are interested only in the creation of the relic gravitational waves by the expanding spacetime 
background, the perturbed matter source is therefore not taken into account. Moreover, as the relic 
gravitational waves are very weak, in the sense that hij <^ 1, so one need just study the linearized field 
equation: 

d^{./^d^h,,{^,T)) = 0. (3) 

In quantum theory of gravitational waves, the field hij is a field operator, which is written as a sum of 
the plane wave Fourier modes 



h„(x, r)^^ ± r <i'k4^)(k)4=[4^'4^)(r)e'-- + hi''- (r)e-"^% (4) 

where Ipi — ^fG is the Planck's length, the two polarizations S^{}^), \ — 1,2, are symmetric and 
transverse-traceless k'ef^i^ = 0, ?>'^ef^{^ = 0, and satisfy the conditions e(^)*^(k)eSf ^(k) = 2(5aa', 
and ej^^(— k) = e-^^(k), the creation and annihilation operators satisfy [ak^\ Ok'^ ^] = ^w^^i^ — k'), 
and the initial vacuum state is defined as 

(^i^ I >= (5) 

for each k and A. As a matter of fact, this definition depends on the choice of the mode function 
h^i^\r); different h^i^\T) define different vacuum states, a point that will be further explained later. 
In the vacuum state the energy density of gravitational waves too = 3^hij,oh^^,o gives the zero-point 
energy 

J 2 J —oo 

For a fixed wave number k and a fixed polarization state A, the Eq.(3) reduces to the second-order 
ordinary differential equation [20] 

A«"(r)+2^ftf'M + ^'A«=0. (6) 

where the prime denotes d/dr. Since the equation of /ik^^(r) for each polarization is the same, we 
denote h^\T) by /ik (''■)■ One rescales the field hi^(r) as 

h.{r)^, (7) 

a(r) 



< 



a" 



and the equation for //k is 

+ ik' - -hk = 0. (8) 

This equation can be regarded as the equation for one-dimensional oscillator in a given effective potential 
barrier a" /a. For a given spacetime background with a generic power-law form of the scale factor 

a(T) (X r", (9) 



where a is a real number, positive or negative, the general solution is a linear combination of the Hankel's 
functions: 



fik{r)AkVk^H[^_,JkT) + BkVk^H[^J_,^{kT). (10) 

Given a model of the expansion of Universe, consisting of a sequence of successive a(r) as in Eq.(9) with 
different a, one can construct an exact solution /ifc(r) by matching its values and its derivatives at the 
joining points. One may may also numerically solve Eq.(6) and give the corresponding power spectrum 
as we will present later in the paper. The vacuum state |0 >, determined by the mode function /^^(t), 
is therefore fixed by a choice of coefficients Ak and Bk. In the limit A; — > oo, or r — > ±00, 



("-2)^ ^ V TT 



approaching the positive and negative frequency modes, respectively. For instance, if we choose = 0, 
then in the limit A; — > 00, or r — > ±00, 

/ik(r) oc e-"'^ 

giving the positive frequency mode. This choice yields the so-called adiabatic vacuum [21] through (4) 
and (5). The amplitude of the power spectrum of the gravitational waves depends only on the initial 
value of |/ik(''')| at the horizon-crossing. 

The following two limiting cases are useful for an approximate evaluation of the spectrum. Outside 
the barrier A;^ ^ ^ (equivalent to k ^) the gravitational wave field reduces to 

/ife(T) ^^fe— - + -, (11) 

a{T) a(r) 

having a decreasing amplitude 

hk(r) oc l/a(T). (12) 
Inside the barrier A;^ ^ ^ (equivalent to k <^ ^) the gravitational wave field 

a^(r'j 

The second term D J'^^,^ is small for the models that we shall study in the following, so the long 
wave-length limit of hk is simply a constant: 

hk{T) = C,. (14) 

Thus, as a function of r, /ifc(r) has simple approximate behaviors in the two limiting cases, and we will 
use these to estimate the spectrum at present stage. 

3. Epochs of Expanding Universe 



The history of overall expansion of the Universe can be modelled as the following sequence of 
successive epochs of power-law expansion. 
The initial stage (inflationary) 

a(T) = Zo I T r+^, -oo < r < (15) 

where 1 + /5 < 0, and ti < 0. The special case of /? = — 2 is the de Sitter expansion. 
The z-stage 

a{T) ^ a,{T - TpY+''% n<T<Ts, (16) 

where /3s + 1 > 0. Towards the end of inflation, during the reheating, the equation of state of energy in 
the Universe can be quite complicated and is rather model-dependent [22]. So this z-stage is introduced 
to allow a general reheating epoch, as has been advocated by Grishchuk [8]. 
The radiation-dominated stage 

a{r) = ae{r - Te), Tg <r < T2. (17) 

The matter-dominated stage 

a{T) = am{r - Tmf, T2<T < Te, (18) 

where te is the time when the dark energy density pA is equal to the matter energy density pm- Before 
the discovery of the accelerating expansion of the Universe, the current expansion was usually taken to 
be in this matter-dominated stage, which is a decelerating one. Now, following the matter-dominated 
stage, we add an epoch of accelerating stage . The value of the redshift ze at the time te is given by 
(1 + Ze) = a{TH)/a{TE), where th is the present time. Since p\ is constant and Pm{T) oc a~^(r), one 
has 

1 = = PA 

Pm{rE) Pm{rH){l + ZeY' 

If we take the current values Q\ ~ 0.7 and ~ 0.3, then it follows that 

1 + ZE ^ 1.33. 

The accelerating stage (up to the present) 

a(r) ^ Ih \ T - Ta \-\ te<t<th. (19) 

This stage describes the accelerating expansion of the Universe, which is the new feature in our model and 
will induce some modifications to the spectrum of the relic gravitational waves. It should be mentioned 
that the actual scale factor function a(T) differs from Eq.(19), since the matter component exists in 
the current Universe. However, the dark energy is dominant, so (19) is an approximation to the current 
expansion behavior. 



Given a(T) for the various epochs, the derivative a' — da/dr and the ratio a' /a follow immediately. 
Except for the Ps that is imposed upon as the model parameter, there are ten constants in the above 
expressions of a(r). By the continuity conditions of a(r) and a(r)' at the four given joining points 
Ti, Ts, T2, and te, one can fix only eight constants. The remaining two constants can be fixed by the 
overall normalization of a and by the observed Hubble constant as the expansion rate. Specifically, we 
put \ Th — Ta \= I as the normalization of a, which fixes the constant r^, and the constant In is fixed 
by the following calculation 

\ / TH 

so Ih is just the Hubble radius at present. Then everything is fixed up. In the expanding Robertson- 
Walker spacetime the physical wave length is related to the comoving wave number by 

A . (21) 
and the wave number kn corresponding to the present Hubble radius is 

kH = = 2n. (22) 



There is another wave number. 



2'iTa{rE) kn 

kR — 



1/H 1 + ze' 

whose corresponding wave length at the time te is the the Hubble radius 1/H. 
By matching a and a' /a at the joint points, we have derived, for example, 

lo = lHbCE^'^^k!^C^CP^, (23) 

where b = \1 + which is defined differently from Grishchuk's [8], (e = C2 = 

(^g = ^1^, and = With these specifications, the functions a(r) and a'{T)/a{T) are fully 

determined. In particular, a'{T)/a{T) rises up during the accelerating stage , instead of decreasing as in 
the matter-dominated stage. This causes the modifications to the spectrum of relic gravitational waves. 

4. The Spectrum of Gravitational Waves. 

The power spectrum h{k,T) of relic gravitational waves is defined by the following equation 



L 



00 Ah 

h\k,T)— =< 0|/i^^(x,T)/ii,-(x,T)|0 >, (24) 

k 



where the right-hand-side is the vacuum expectation value of the operator Whij. Substituting Eq.(4) 
into the above, and taking the contribution from each polarization to be the same, one reads the power 
spectrum 

h{k,T)^^k\hk{T)\ (25) 



Once the mode function /ife(r) is given, the spectrum h{k,T) follows. 

The initial condition is taken to be during the inflationary stage. For a given wave number k, its 
wave crossed over the horizon at a time Tj, i.e. when the wave length Aj = 27ra{Ti)/k is equal to 
l/H{Ti), the Hubble radius at time n. Eq.(15) yields Hiji) = lo\Ti\'^^^ /\1 + f3\, and for the exact de 
Sitter expansion with P = —2 one has H(Ti) = Iq. Note that a different k corresponds to a different 
time Tj. Now choose the initial condition of the mode function hkij) as 

\hu{n)\ = ^ (26) 



Then the initial amplitude of the power spectrum is 



M^,T,) = 8v^^. (27) 
From Aj = l/H{Ti) it follows that = ^. So the initial amplitude of the power spectrum is 

h{Kn)^A{^f+^ (28) 

A = 8^/^6^. (29) 
'o 

For the case of /3 = — 2 the initial spectrum is independent of k. The power spectrum for the primordial 

perturbations of energy density is P{k) oc \h{k, th)\'^, and its spectral index n is defined as P{k) oc /c"~^. 

Thus one reads off the relation n = 2/? + 5. The exact de Sitter expansion with (5 — —2 will yield the 

so-called scale-invariant spectral index n — 1. 

Once the initial spectrum is specified, we can derive the spectrum h{k,TH) at the present time th- 
For k < kE, the wave lengths in this range are even greater than the present Hubble radius 

Ih, one has k < a'/a throughout the whole expansion up to the present, so the amplitude remains the 

same constant as the initial one in Eq.(28): 

h{k,TH)=A{-l^f+^, k<kE. (30) 

kn 

During the whole period inside the barrier, the spectral amplitude remains h{k,Ti) approximately 
until the wave leaves the barrier and begins to decrease as l/a(r). Let a^^{k) be the scale factor at 
this moment. For those very long wavelength modes with kE < k < kn, during the current epoch of 
accelerating expansion, /ifc(r) stops decreasing as soon as the barrier a'/a is higher than k at a time r(/c) 
earlier than th, so hk{T) has decreased by a factor and the amplitude of the present spectrum 

is given by 

KKrH) = A{^r^^f^ (31) 
ku a{T{k)) 

The decreasing factor is written as 

a**(fc) a^^{k) aijE) 



where the constant 



a{T{k)) a(TE) a(r(k))' 



During the matter-dominated stage one has a**(/c) oc oc A; ^ and a(rE) oc rl oc k^^, and during the 
accelerating stage one has a{TE) oc r^^ oc and a(r(A;)) oc oc k, so the decreasing factor is 



a{T{k)) k 

Thus one gets 

(32) 

where the relation |^ = j^pj^ has been used. The spectrum has a rather stiff slope with the power-law 
index P — 1. The occurrence of this segment of power spectrum is a new feature of the model of 
accelerating expansion that is absent in the decelerating model. The wave lengths corresponding to this 
{kE, kff) segment are very long, comparable to the present Hubble radius, and can only possibly be 
observed through the CMB anisotropies at low multipoles. 

For all the wave number k > kn, as soon as the waves leave the barrier at a^^{k), the modes /ife(r) 
have been decreasing ail the way up to the present time th- So it has been reduced by a factor 
and the amplitude of present spectrum is given by 



h{k,rH)^A{^rP^. (33) 
kn a{TH) 

We use this formula to obtain the following result for the spectrum of the gravitational waves in the 
remaining range of wave numbers. 

For kn < k < k2, the wave number does not hit the barrier, so 

a^*{k) a**(/c) a{TE) _ ^^^2 1 



airn) a{rE) a{rH) k (1 + ze) ' 
one obtains 

h{k,rH) ^ A{-^f ^ kH<k<k2. (34) 

Kh U + Ze)'^ 

The spectrum in this interval differs from that of the matter-dominated model by an extra factor 
(^l_^.],^)3 — ^ QA3. The wave-lengths in this range are very long, but are still shorter than Ih- The 
spectrum in this interval may contribute to, and, therefore, have its imprints in CMB anisotropies. Let 
us estimate the value k2. Assuming that the equality of radiation and matter occurred at the redshift 
Z2 ^ 3454, as indicated by the WMAP observation [19], one has 

h ^ .a(r^a ^ Mth} , a{r^ . ^ ^3454^ ~ 51 
For k2 < k < kg, the calculation is similar to the previous case with the result 

h{k,TH) = Mh'^'i^)7T^^,, k2<k<K. (35) 

kn k2 (1 + zeY 



Again the extra factor 1/(1 + appears. Notice that this range of frequency covers the one that the 
detectors of LIGO and LISA are sensitive. 
For kg < k < ki, the calculation yields 

KKm) = K<k<k,. (36) 

kn kn k2 (1 + zsr 

It also contains an extra factor (jq:^^- This is the high frequency range which is still beyond the current 
detection. 

5. Determining the Parameters. 

To completely determine the spectrum, we also need to specify the values of ui, A, Ug, Ps, that appear 
in the expressions for the spectrum h{k,TH)- Since the wave number is proportional to the frequency, 
k cc u, the ratios of the wave numbers can be replace by that of the frequencies, e.g., k/ku = v jvu, etc., 
in the above formulae of the spectrum h{k^ th)- The Hubble frequency uh = 1/Ih = H 2x 10~^*Hz. 

An estimate of the highest frequency ui can be made as given in [8]. From the expression (36) one 

has 

M^i,r^) = 8V^^(^) = 8V^^, (37) 

where the expressions (23) for b/lo and (29) for A have been used. The spectral energy density parameter 
ilg{iy) of the gravitational waves is defined through the relation pg/pc = JO,g{u)^, where pg is the 
energy density of the gravitational waves and pc is the critical energy density. One reads 

n,H-'^h\k,rH)(^r- 

6 Uh 

If it is imposed that at the highest frequency ui the value ^g{i^i) not exceed the level of 10~^, as required 
by the rate of the primordial nucleogenesis, then one gets ui — 3 x lO^^Hz. 

Next let us estimate the overall factor A in the spectrum h{k,TH)- If the CMB anisotropies at low 
multipoles are induced by the gravitational waves, or, if the contributions from the gravitational waves and 
from the density perturbations are of the same order of magnitude, we may assume AT/T ~ h{k, th)- 
This will determine A. The observed CMB anisotropies [19] at lower multipoles is AT/T ~ 0.37 x 10~^ 
at Z ~ 2, which corresponds to the largest scale anisotropies that have observed so far. Taking this to 
be the perturbations at the Hubble radius 1/H yields 

h{kH, Th) = A-—^—-r = 0.37 x 10"^ (38) 

Actually the observed AT/T varies with I, for instance, AT/T ~ 1.04 x 10~^ at I ~ 10, which would 
give an outcome for A greater than that in Eq.(38) by a factor ~ 2.8. 

However, there is a subtlety here on the interpretation of AT/T at low multipoles, whose corre- 
sponding scale is very large ~ Ih- At present the Hubble radius is Ih, and the Hubble diameter is 21h- 



On the other hand, the smallest characteristic wave number is kE, whose corresponding physical wave 
length at present is ^""^^^^^ = ///(I + ze) — I.S21h, which is within the Hubble diameter 21h, and 
is theoretically observable. So, instead of (38), if AT/T ~ 0.37 x 10~^ at Z ~ 2 were taken as the 
amplitude of spectrum at ue, one would have h{kE,TH) = ^(jqi^^ji+F — 0-37 x 10~^, yielding a smaller 
A than that in Eq.(38) by a factor (1 + ^e)^"^ ~ 2.2. 

We now check the range of P allowed. During the inflationary expansion when the k-mode wave 
entered the barrier with Aj = l/Hiji), from which it follows that Xi = ^(^)^^^. For the classical 
treatment of the background gravitational field to be valid, this wave length should be greater than the 
Planck length, > Ipi, so 

A ■ 

At the highest frequency v — vi, this yields a constraint 

/3<-2 + ln(^)/ln(^), (39) 

which depends on A. Thus, for A given in (38), one obtains the upper limit (3 < —1.78. We remark 
that this range is larger than that in the decelerating model [8], which allows for only P < —1.9. 

Finally we give an estimate of the allowed values of Ps and Ug. Plugging b/lo given by (23) into 
A = using U2/uh = 58.8 and Ih/Ipi = 1-238 x 10^^, one has 



1.484 X 10^ ^ ^ = (i^)-/3(^)/3.. (40) 



, 7^ = (^)-^(^- 

,1 + Ze)-^ fH l^s 

Given a set values of A and /?, one can take Vs and Ps to satisfy this relation. From Eq.(36) it is seen 
that, for a fixed z/^, a smaller Ps tends to increase slightly the amplitude hiy^TH)- For definiteness, we 
take Vg = 10^ Hz as an example in the following. For A given in (38), taking P — —1.8, —1.9, then 
ps = 0.598, -0.552, respectively 

With these results we plot the spectrum h{k,TH) vs the frequency v for two values of parameter 
P — —1.9, —1.8 in Figure 1. Besides, for convenience of comparison, the spectrum for the decelerating 
model is also plotted for P — —1.9. It is clearly seen that, for the same parameter P — —1.9 in the 
whole range v > uh, the accelerating model has an amplitude lower by a factor of ~ 2.2 than that 
of the decelerating model. This is due to the extra factor ^^^^^^^^3 = as has been demonstrated in 
the previous expressions (34), (35), and (36). Interestingly, the larger value P — —1.8 gives a flatter 
spectrum h{i',TH) with an overall higher amplitude in the range v > 1/2. In particular, in the higher 
frequency range (10^^ ~ 0) Hz covered by LISA, the amplitude is /i ~ 10~^^ ~ 10~^^ for the P — —1.8 
case, about 10 to 10^ higher than the P — —1.9 case. In the even higher frequency range (10 ~ 10^) 
Hz covered by LIGO, the amplitude is /i ~ lO^^^ ~ 10"^^ for the P = -1.8 case, also 10 to 10^ higher 
than the P — —1.9 case. Thus an inflationary model of larger index n — 2p + 5 will predict a stronger 
signal of relic gravitational waves in higher frequencies. 
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Figure 1: The spectrum h{k, Th) as a function of the frequency v in the present Universe in both 
accelerating and decelerating expansion. 
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Figure 3: The energy spectrum ^^(z/) in the present Universe in both accelerating and decelerating 
expansion. 

The recent second science run of the LIGO interferometric detectors [23] gives a sensitivity 3 x 10"^'' 
to 3 X 10^^"^ in the frequency range (10^, 10'^)Hz. The best sensitivity is given by the 4km arm LI 
detector located at Livingstone, which is about 3 x 10"^^ near a frequency ~ 300Hz. Our calculation for 
the P = —1.9 case yields an amplitude, h ~ 10^^^, much smaller than this sensitivity of the second run. 
However, the most interesting case is the model of P = —1.8, in which the amplitude of the gravitational 
waves is just to fall into the sensitivity of the LI detector. Since the second run of LIGO has not observed 
any signal of stochastic gravitational waves in this frequency range, we arrive at a constraint P < —1.8 
on the model parameter of the inflationary expansion. Note that this constraint is also consistent with 
P < —1.78, as has been imposed, through Eq.(38), from the observed CMB anisotropics of WMAP. 

As a double check, we have also numerically solved the differential equation (6) of the gravitational 
waves, and have found the resulting power spectrum h(u, t) from the numerical solution, which is plotted 
in Figure 2. Since the numerical result that we have plotted carries the oscillating factor cos(27rz/(t — t,^)), 
so its curve shows an extra small zigzag, as expected. Moreover, for the range of small frequencies, 
u < uh, cos(27rz/(t — t^)) ~ 1, the oscillating amplitude is very small, confirming an observation made 
by ref. [8]. Except for this oscillating effect, the numerical result of the spectrum in Figure 2 agree with 
the analytical one in Figure 1. 

We also plotted the spectral energy density from the analytic solution in Figure 3. It is 

similar to the known result, except for the obvious distortions caused by the acceleration of the Universe 
expansion in the low frequency range . 




6. Conclusion. 



We have presented a calculation of the spectrum of relic gravitational waves in the present Universe in 
accelerating expansion. The recent WMAP result of AT/T has been used to normalize the amplitude of 
relic gravitational waves. In comparison with the decelerating models, the spectrum has been modified 
due to the particular form of the barrier function a' /a during the acceleration of current expansion. 
Specifically, in very low frequency range (< vh) the spectrum has been changed with the peak of 
spectrum being now located at ve, and there appears a new segment of spectrum from ve to vh- 
These very long wave length features can be only be possibly detected by the CMB anisotropies at low 
multipoles. In the higher frequency range (> vh) the spectral amplitude acquires an overall factor ^ 
as compared with the decelerating model. This higher frequency range is pertinent to the detection 
projects such as LISA and LIGO. In addition, using the AT/T at low multipoles for normalization gives 
a larger upper limit P < —1.78. A larger value of (3 yields a flatter spectrum h{v,TH) as a function of 
V, producing more power on the higher frequencies. The resulting sensitivity of the second scientific run 
of the LIGO detectors, has put a restriction on the model parameter (3 < —1.8. 
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